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Abstract: Determining the mechanisms controlling the changes of wet and dry conditions will improve 
our understanding of climate change over the past hundred years, which is of great significance to the 
study of climate and environmental changes in the arid regions of Central Asia. Forest trees are 
ecologically significant in the local environment, and therefore the tree ring analysis can provide a clear 
record of regional historical climate. This study analyzed the correlation between the tree ring width 
chronology of Juniperus turkestanica Komarov and the standardized precipitation evapotranspiration 
index (SPEI) in Northwest Tajikistan, based on 56 tree ring samples collected from Shahristan in the 
Pamir region. Climate data including precipitation, temperature and the SPEI were downloaded from the 
Climate Research Unit (CRU) TS 4.00. The COFECHA program was used for cross-dating, and the 
ARSTAN program was used to remove the growth trend of the tree itself and the influence of non- 
climatic factors on the growth of the trees. A significant correlation was found between the radial growth 
of J. turkestanica trees and the monthly mean SPEI of February—April. The monthly mean SPEI 
sequence of February—April during the period of 1895-2016 was reconstructed, and the reconstruction 
equation explained 42.5% of the variance. During the past 122 a (1895-2016), the study area has 
experienced three wetter periods (precipitation above average): 1901-1919, 1945-1983 and 1995-2010, 
and four drier periods (precipitation below average): 1895—1900, 1920-1944, 1984-1994 and 2011- 
2016. The spatial correlation analysis revealed that the monthly mean SPEI reconstruction sequence of 
February—April could be used to characterize the large-scale dry-wet variations in Northwest Tajikistan 
during the period of 1895-2016. This study could provide comparative data for validating the 
projections of climate models and scientific basis for managing water resources in Tajikistan in the 
context of climate change. 
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1 Introduction 


Climate change is a global scientific problem that is directly related to human survival, the 
environment and the formation of an ecological civilization (Immerzeel, 2010). There are 
significant regional differences in climate change, with changes in precipitation being especially 
local. Thus, it is particularly important to understand the consistency of regional dry and wet 
periods. Central Asia is a typical inland arid region that is located in the confluence of large-scale 
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atmospheric circulation systems that have shown strong responses to global climate change (Wolff 
et al., 2017). The environment in Central Asia is fragile and is sensitive to climate change. 
Therefore, it is essential to understand the climate history of this region when studying global 
climate change and the sustainable use of regional water resources (Seddon et al., 2016). 

Tree rings have been widely used to reconstruct historical climate and their use has been shown 
to have high resolution, replicability, strong continuity, a wide sample distribution and accurate 
dating (Briffa et al., 1998; Cook et al., 2010; Liu et al., 2010; Hughes et al., 2011; Pederson et al., 
2014; Gou et al., 2015; Zhang, 2015; Wang et al., 2016; Zhang et al., 2018). The standardized 
precipitation evapotranspiration index (SPEI) is calculated by introducing potential evaporation 
factors to the standardized precipitation index (Vicente-Serrano, 2010). The SPEI is a relatively 
new drought index with the advantage of multiple time scales and considering the influence of 
temperature changes on drought (Wang and Chen, 2014). The SPEI can not only identify the 
beginning and end time of a drought event, but also measure its severity according to drought 
intensity and duration. The index has been widely used for drought analysis and assessment 
worldwide (Allen et al., 2011; Vicente-Serrano et al., 2011; Bao et al., 2015). 

Although there are many climate reconstruction sequences in Central Asia, most are concentrated 
in the Tianshan Mountains. There are few tree ring sites in the Pamir region in Tajikistan and 
therefore few climate reconstruction sequences have been reported. Opata-owczarek et al. (2014) 
proposed that shrubs from the Pamir region in Tajikistan have the potential for use in tree ring 
climatology research, and then they produced a 1070-a climate time series through the crossover of 
dead and live trees of the same species in the Pamir region (Opata-owczarek et al., 2018). Opata- 
owczarek et al. (2019) also suggested that the Pamir region in Tajikistan is an ideal area to 
reconstruct historical climate during the past millennium. 

Here, we presented a tree ring based drought variability in the Pamir region of Northwest 
Tajikistan, spanning the period from 1895 to 2016. This study reconstructed the SPEI with tree ring 
data, which were used to identify dry and wet periods in Northwest Tajikistan over the past 122 a. 
The historical monthly mean SPEI of February—April during the period of 1895-2016 were 
reconstructed. Based on the reconstruction results, we analyzed the characteristics of historical dry- 
wet variations and revealed the main trends of dry and wet conditions in this area. This study could 
provide comparative data for validating the projections of climate models in Central Asia. In 
addition, the results of the study have great significance for understanding the evolution of the 
climatic environment in Central Asia under the background of global climate change. 


2 Materials and methods 


2.1 Study area 


Tree ring samples were collected from Shahristan (39°34'N, 68°50'E; 2600-2800 m a.s.l.) in the 
Pamir region of Northwest Tajikistan. The Pamir region spans Tajikistan, China and Afghanistan 
(Qiu, 2005), and is the confluence of many major mountain ranges in Central Asia. The study area 
is characterized by a temperate climate, with average annual precipitation of 500 mm. The monthly 
precipitation is largest in spring (March, April and May). Summer is a period of high temperatures, 
with the highest average temperature of 16.0°C in July. The lowest average temperature is —11.8°C 
in January. Juniperus turkestanica Komarov, a dominant evergreen species with strong adaptability 
in this region, is highly tolerant to cold and drought. This species can grow on mountain ridges or 
riparian areas. 


2.2 Climate data collection 


Precipitation, temperature and SPEI data in the study area were downloaded from the Climate 
Research Unit (CRU) TS 4.00 (http://www.cru.uea.ac.uk) (Harris et al., 2014), with a spatial 
resolution of 0.5°x0.5° for the time period of 1930-2016. The SPEI has been widely used for 
monitoring the spatial and temporal variability of drought and its changing tendency in many 
different regions. To determine the correlation between the tree ring width chronology and climate 
factors, we used CRU data with a spatial resolution of 0.5°x0.5° in the spatial correlation analysis. 
Specifically, mean precipitation and temperature data of grid points (38°75'-41°75’N, 66°75'— 
72°25'E) during the period of 1930-2016 were used in the spatial correlation analysis. The 
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downloaded data were processed by using MATLAB. 
2.3 Tree ring sample collection and tree ring width data 


In October 2016, 60 tree ring samples were collected from 30 healthy J. turkestanica trees in 
Shahristan in the Pamir region of Northwest Tajikistan. The height of the trees sampled was 
between 6.0 and 13.0 m, and the diameter was between 80.0 and 210.0 cm. Generally, sampling 
was conducted at breast height (about 1.3 m) and two to three 5.0 mm diameter cores from one tree 
were extracted with borers applied from various directions for further dendrochronological dating. 
However, the crown of some trees was too low for sampling at breast height; thus, samples were 
collected at a height of 0.2-0.5 m above the ground. The cores were removed and placed in 
cardboard tubes. The inner diameter of the cardboard tube was 1.0 mm wider than that of the tree 
ring sample, and once the sample was in place, it was sealed with tape and the cardboard was 
labelled. 

The samples were brought back to the laboratory in strict accordance with established tree ring 
analysis procedures (Fritts, 1976). The samples were dried, fixed and polished until the spring 
wood and summer wood could clearly be distinguished under the microscope. Then, cross-dating 
was conducted and the width of the tree rings was measured. A tree ring measuring instrument 
(0.01 mm) was used to measure the width of the tree rings, and the COFECHA program was used 
for cross-dating (Holmes, 1983). It should be noted that we removed the samples that were weakly 
correlated with the main sequence according to the COFECHA results. Due to dating errors 
connected with missing rings, wedging rings, false rings and other wood anomalies, four samples 
were removed. Finally, a total of 56 samples from 28 healthy J. turkestanica trees were used for 
building the tree ring width chronology. The ARSTAN program was used to remove the growth 
trend of the tree itself and the influence of non-climatic factors on the growth of the trees (Cook, 
1985). The double weighted average method was used to standardize the tree ring sequence. Then, 
the autoregressive and standardized tree ring width chronologies were obtained. The standardized 
tree ring width chronology was further used for the correlation analysis. 

To identify the main limiting factors for the growth of tree ring width in the study area, we used 
DendroClim2002 (Biondi and Waikul, 2004) to determine the strength of the correlation between 
the tree ring width chronology and climate data from June of the previous year to September of the 
current year. Then, we established a regression equation based on the results of the correlation 
analysis between the tree ring width chronology and climate factors. We also applied the 
elimination test method to evaluate the effectiveness of climate reconstruction. The correlation 
coefficient and standard deviation were calculated, and the variance interpretation, a sign test and 
the reduction of error test were performed (Cook and Kairiukstis, 1990). 


3 Results 
3.1 Tree ring width chronology 


Because the standardized tree ring width chronology retains more low frequency information than 
the autoregressive tree ring width chronology, it can be more responsive to natural climate change. 
The standardized tree ring width chronology was therefore used for the analyses (Fig. 1). The 
sequence length of the standardized tree ring width chronology was 316 a (1701-2016), and the 
length of the reliable chronological sequence (expressed population signal>0.850) was 122 a 
(1895-2016). The mean sensitivity was 0.263, the standard deviation was 0.291, and the variance 
explained by the first principal component was 42.5% (Table 1). This indicates that the tree growth 
in the study area was sensitive to climate change. The signal-to-noise ratio was 5.946, and the 
expressed population signal was 0.856, which exceeded the minimum threshold of 0.850 (Efron, 
1979). These results showed that the tree ring width chronology could be considered reliable data 
series for studying climate responses. 
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Fig. 1 Standardized tree ring width chronology and the amount of tree ring samples (i.e., sample depth) 


Table 1 Statistical characteristics of the standardized tree ring width chronology 


Statistic Standardized tree ring width chronology 
Number of sampled trees 28 
Number of tree cores 56 
Chronology period 1701-2016 (316 a) 

R 0.373 

AGR 2.180 

MS 0.263 

SD 0.291 

ACI 0.498 

SNR 5.946 

VFE 0.425 

EPS 0.856 

Year with EPS>0.850 1895 


Note: R, all-series Rbar; AGR, average growth rate; MS, mean sensitivity; SD, standard deviation; AC1, first autocorrelation; SNR, signal- 
to-noise ratio; VFE, variance explained by the first principal component; EPS, expressed population signal. 


3.2 Relationship between the radial growth of tree rings and climate factors 


Considering the lag effect of climate factors on tree growth, June of the previous year to September 
of the current year were selected as the relevant months for analysis. Correlation analyses were 
conducted between the tree ring width chronology of J. turkestanica and the monthly precipitation, 
monthly mean temperature and monthly SPEI (Fig. 2). The result showed that the radial growth of 
J. turkestanica was positively correlated with precipitation and temperature. The radial growth of J. 
turkestanica was significantly positively correlated with the precipitation in November of the 
previous year and in April, July and August of the current year (P<0.05). Moreover, the radial 
growth of J. turkestanica was significantly positively correlated with the temperature in December 
of the previous year and in January of the current year (P<0.05). It was also found that the 
correlation between the radial growth of J. turkestanica and seasonal combination of precipitation 
and temperature was generally low. Therefore, the SPEI was used in the correlation analysis in this 
study. The results showed that the correlation between the radial growth of J. turkestanica and the 
monthly mean SPEI of February—April in the current year was maximal, with the correlation 
coefficient of 0.42 (P<0.01). 


3.3 Reconstructed SPEI of 1895-2016 


The standardized tree ring width chronology was an independent variable. The monthly mean SPEI 
of February—April was downloaded from the CRU TS 4.00 (spatial resolution of 0.5°x0.5°) and 
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Fig. 2 Correlation coefficients between the tree ring width chronology of J. turkestanica and monthly climate data 
during the period of 1930-2016. SPEI, standardized precipitation evapotranspiration index. For the x-axis, p6—p12 
denote June to December of the previous year, respectively; and 1—9 denote January to September of the current 
year, respectively. The dotted lines denote 99% and 95% significance levels. 


used as the dependent variable. The reconstruction equation was established by a linear regression 
model: 
y=1.181x+1.540, (1) 
where y is the monthly mean SPEI value of February—April and x is the standardized tree ring 
width chronology. This reconstruction equation could explain 42.5% of the monthly mean SPEI 
(February, March and April) chronology variance, and the variance interpretation after adjusting for 
degrees of freedom was 40.5%. The correlation coefficient was 0.425 (P<0.01). The synchronicity 
between the reconstructed and downloaded values was sufficient, and therefore the reconstruction 
equation could be used to reconstruct the monthly mean SPEI of February—April in the current year. 
The stability of the equation directly affected the quality of the reconstructed sequence. The 
elimination method was used to cross-check the reconstruction equations and determine their 
stability and reliability. The test statistics are shown in Table 2. The sign test was significant 
(P<0.01), indicating that the reconstructed SPEI sequence was consistent with the downloaded 
sequence in terms of the high and low frequency variations. The product mean test value was 3.73 
(P<0.01). The reduction of error value was 0.15 (positive value), which means that the climate 
estimates passed the test. Therefore, the equation used to reconstruct the SPEI in Northwest 
Tajikistan is reliable. 


Table 2 Statistical characteristics of a leave-one-out test 
r r SD RE/CE t ST F 


0.425 0.181 49.12 0.15/0.15 3.73 57+/29- 16.31 


Note: r, correlation coefficient; 7°, test inter-series correlation; SD, standard deviation; RE, reduction of error; CE, coefficient of efficiency; 
t, product mean test; ST, sign test; F, regression equation test. 


The SPEI sequence of February—April during the period of 1895-2016 was reconstructed, and 
the result is shown in Figure 3. An 11-a moving average of the SPEI was also reconstructed (Fig. 4) 
to determine the characteristics of inter-decadal dry-wet variations in the study area. The study area 
experienced four dry periods (precipitation below average): 1895—1900, 1920-1944, 1984-1994 
and 2011-2016 during the past 122 a. The driest period was 1920-1944, with the monthly mean 
SPEI value being 0.24 below the average of 1895-2016. The longest dry period lasted for 25 a 
(1920-1944; Table 3). During the past 122 a, the study area also experienced three wet periods 
(precipitation above average): 1901-1919, 1945-1983 and 1995-2010. The wettest period was 
1901-1919, with the monthly mean SPEI value being 0.52 above the average of 1895-2016. The 
longest wet period lasted for 39 a (1945-1983; Table 3). During the study period, the driest year 
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was 1937 and the wettest year was 1907. The results of the SPEI reconstruction showed that there 
has been a drying trend in Northwest Tajikistan over the last six years (2011-2016; Table 3). 
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Fig. 3 Comparison of the downloaded and reconstructed monthly mean SPEI of February—April mean SPEI for the 
period of 1930-2016 
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Fig. 4 Reconstructed monthly mean SPEI of February—April and the 11-a moving average for the period of 1895- 
2016 


Table 3 Wet-dry periods based on the standardized precipitation evapotranspiration index (SPEI) values for the 
period of 1985-2016 


Dry period Number of Anomaly Wet period Number of Anomaly 
years (a) years (a) 

1895-1900 6 0.09 1901-1919 19 0.52 

1920-1944 25 0.24 1945-1983 39 0.11 

1984-1994 11 0.08 1995-2010 16 0.18 

2011-2016 6 0.20 


Note: Anomaly was calculated by the mean SPEI of this period divided by the mean SPEI of the study period (1895-2016). 
3.4 Spatial representation of SPEI 


To determine the relationship between the SPEI reconstruction and the geographical representation 
(Fig. 5), we analyzed the spatial correlation between the reconstructed monthly mean SPEI of 
February—April and the downloaded monthly mean SPEI of February—April from the CRU TS 4.00 
(spatial resolution of 0.5°x0.5°) for the period of 1930-2016. The results of this reconstruction for 
1930-2016 suggested broad-scale regional climatic variations. The spatial correlation analysis 
revealed that the reconstructed SPEI sequence could depict the large-scale dry-wet variations in 
Northwest Tajikistan. 
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Fig. 5 Spatial correlation fields of (a) reconstructed monthly mean SPEI of February—April and (b) downloaded 
monthly mean SPEI of February—April from the Climate Research Unit (CRU) TS 4.00 in Northwest Tajikistan for 
the period of 1930-2016 


4 Discussion 


4.1 Radial growth of J. turkestanica response to climate 


Numerous studies have shown that the radial growth of trees is not only influenced by the climate in 
the growing season, but also by the seasonal climate before the growing season (Fritts, 1976). The 
radial growth of trees is controlled by the moisture conditions in the early growing season (Yuan et 
al., 2001, 2003; Zhang et al., 2016). The radial growth of trees benefits from precipitation during the 
early growing period because it enhances soil water storage, which is important for tree growth in 
the following year. Precipitation in the growing season directly increases soil water availability, thus 
compensating for the evapotranspiration of soil water (Li et al., 2006). In the early growth stage of J. 
turkestanica trees, the moisture required for tree growth is mainly dependent on the melting of snow 
cover, and therefore snowmelt water in spring has an important influence on the radial growth of 
trees. At the monthly scale, February to April is a crucial period in which trees enter the growing 
period from dormancy and moderate warming could help to break their dormancy early. High 
temperatures during February—April can promote snow melt, increasing the soil moisture content 
and promoting the radial growth rate of trees in the growing season (Li et al., 2011). 

The effects of temperature changes on drought have been evaluated using the SPEI (Wang and 
Chen, 2014). The SPEI can not only identify the beginning and ending time of a drought event, but 
also measure drought severity according to its intensity and duration. In Northwest Tajikistan, there 
was a positive relationship between the radial growth of J. turkestanica trees and the SPEI, as 
shown in Figure 2. This finding is consistent with the result in eastern Xinjiang of China (Chen et 
al., 2015). The larger the SPEI, the higher the degree of wetness; in contrast, the smaller the SPEI, 
the stronger the degree of drought. The correlation coefficient between the standardized tree ring 
width chronology and the monthly mean SPEI of February—April was high (7=0.425, P<0.01). 
February to April is the early growing season of J. turkestanica, which is the key period for the 
early wood formation of trees. Wet conditions in this period will therefore promote the rapid 
growth of early wood (Chen et al., 2011). 


4.2 Implication of the reconstructed SPEI 


The study area has experienced significant humidification throughout the 20" century. During the 
study period, a total of 64 a that was considered wetter than the average and a total of 37 a that was 
drier than the average. Reconstructions of precipitation in Northwest China and evaluations based 
on the Palmer drought severity index (PDSI) have indicated a significant drying trend in recent 
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decades (Li et al., 2007; Fang et al., 2012; Liu et al., 2013a, b; Yang et al., 2016). These results are 
consistent with the finding of our study. In Northwest Tajikistan, there was a trend toward drought 
conditions after 2011, which was especially evident in the recent six drier years (2011-2016). 
Based on a tree ring chronology reconstruction in eastern Xinjiang and the northeastern to northern 
part of Tibetan Plateau in China, Liang et al. (2006) found a severe drought occurred at the end of 
the 1920s. Evidence from the reconstruction of precipitation changes at Taibai Mountain in the 
Qinling Mountains and the reconstruction of temporal and spatial precipitation changes in the 
middle and eastern part of Northwest China over the past 400 a showed that the drought was most 
obvious from the 1920s to 1930s (Li et al., 2012; Liu, 2016). These results are consistent with those 
presented in Table 3 that were observed in this study. The reconstructed sequences were spatially 
correlated with the downloaded SPEI datasets, which indicated that the reconstructed SPEI 
sequence had a strong spatial representation. The dry-wet change sequence can better represent the 
dry-wet variations in most areas of Tajikistan. 


5 Conclusions 


Juniperus turkestanica tree ring samples from Northwest Tajikistan were used to build 
autoregressive and standardized tree ring width chronologies. The correlations between the 
standardized tree ring width and the monthly precipitation, monthly mean temperature and monthly 
SPEI were analyzed. The radial growth of J. turkestanica was positively correlated with 
precipitation, temperature and the SPEI. The strongest correlation was found for the relationship 
between the radial growth of J. turkestanica and the monthly mean SPEI of February—April in the 
current year. The study area was found to have experienced three wetter periods (1901-1919, 
1945-1983 and 1995-2010) and four drier periods (1895—1900, 1920-1944, 1984-1994 and 2011- 
2016). A spatial correlation analysis revealed that the reconstructed SPEI sequence could 
characterize the large-scale dry-wet variations in Tajikistan. 

The tree ring width chronology is widely used to reconstruct the regional historical climate, 
however, the stable carbon isotope ratio in late cellulose has been proved one of the most reliable 
indicators of the paleoclimate, with different alternative materials having different responses to 
climate. Therefore, assessments using other tree materials rather than tree rings are necessary to 
accurately reconstruct the historical climate. 
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